import numpy as np

def latlon_to_local(lat, lon, ref_lat, ref_lon):
    R = 6371000  # 地球平均半径（米）
    dlat = np.radians(lat - ref_lat)
    dlon = np.radians(lon - ref_lon)

    # 简化投影（小范围内适用）
    x = R * dlon * np.cos(np.radians(ref_lat))
    y = R * dlat
    return x, y


def local_to_latlon(x, y, ref_lat, ref_lon):
    R = 6371000  # 地球平均半径（米）

    # 反解坐标变换
    dlat = y / R
    dlon = x / (R * np.cos(np.radians(ref_lat)))

    lat = ref_lat + np.degrees(dlat)
    lon = ref_lon + np.degrees(dlon)
    return lon,lat